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ABSTRACT 


TEMPERATURE  CONTROL  OF  A CHEMICAL  REACTOR 


This  paper  presents  the  results  of  applying  various 
controllar  tuning  methods  to  the  feedback  tempera- 
ture control  of  a simulated  chemical  reactor. 
Specifically  the  minimum  error-integral  correlations 
presented  by  Lopez  and  by  Rovlra  are  compared  with 
the  Zlegler-Nlchols  method  and  the  loop-compensa- 
tion method  previously  presented  by  these  authors. 
Based  on  speed  of  response  and  control  stability 
the  best  all-around  techniques  are  shown  to  be 
Rovira's  set-point  tuning  and  the  loop-compensation 
method.  The  closed-loop  responses  of  the  reactor 
temperature  to  changes  In  set-point  and  load  are 
presented. 

INTRODUCTION 

The  modem  control  engineer  finds  that  his  task  of 
designing  feedback  control  systems  Involves  the 
specification  of  the  variables  to  be  measured,  the 
streams  to  be  manipulated,  the  sensors  and  trans- 
mitters to  do  the  measuring,  the  control  valve, 
complete  with  actuator  and  positioner  If  so 
required,  and  the  feedback  controller.  A key  step 
in  the  complete  operation  is  the  tuning  of  the 
controller  parameters  to  the  dynamic  characteris- 
tics of  the  process.  The  success  or  failure  of 
the  design  process  may  depend  on  It. 

The  process  control  literature  of  the  last  four 
decades  contains  a number  of  tuning  methods  and 
correlation  formulas  for  feedback  controllers.  Of 
these,  we  have  chosen  to  consider  In  this  paper 
the  minimum  error  Integral  correlations  of  Rovlra 
[1]  for  set-point  Inputs  and  of  Lopez  [2]  for  dis- 
turbance Inputs,  the  pioneer  method  of  Ziegler  and 
Nichols  [3]  and  the  loop-conqiensatlon  method  pro- 
posed by  Martin  et  al  [A].  These  methods  are 
chosen  because  they  are  Intended  to  apply  In  gen- 
eral to  the  most  comann  processes.  Other  methods 
have  been  presented  In  the  literature  that  apply 
to  more  specific  processes  such  as  the  control  of 
pH,  centrifugal  compressors,  liquid  level,  etc. 

The  methods  studied  here  would  in  no  way  compete 
with  the  more  specific  methods  and  correlations 
for  these  process  Voops. 
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In  order  to  study  the  performance  of  the  com  toller 
to  a more  realistic  plant  than  a first  or  second- 
order  differential  equation,  the  temperature  con- 
trol of  a continuous  stirred  tank  chemical  reactor 
was  chosen.  The  reactor,  which  Is  sketched  in 
Figure  1,  was  simulated  on  an  analog  computer  with 
enough  detail  as  to  Include  Its  most  Important  non- 
linear dynamic  characteristics.  A detailed  descrip- 
tion of  the  mathematical  model  has  been  given  in  a 
previous  publication  [5]. 

The  temperature  In  the  reactor  Is  controlled  *i 
the  desired  set-point  by  manipulation  of  the  late 
of  cooling  water  to  the  jacket.  Step  changes  In 
temperature  set-point  snd  load  are  applied  to 
observe  the  response  of  the  controller  tuned  by 
the  methods  under  consideration.  The  load  Input 
Is  the  rate  of  reactant  feed. 

The  nonlinear  behavior  of  the  loop  Is  lllust  rated 
In  Figure  2 by  a steady-state  plot  of  temperature  - 
the  control  variable  - versus  cooling  water  rate  - 
the  manipulated  variable.  The  slope  of  this  line 
Is  proportional  to  the  process  gain  and  la  shown 
to  Increase  as  the  cooling  water  rata  decreases. 

MODELS  FOR  CONTROLLER  TUNING 

All  of  the  tuning  methods  considered  here  make  use 
of  correlations  which  are  based  on  the  parameters 
of  either  a first-order  lag  plus  dead-time  (trans- 
portation lag  or  time  delay)  or  a second-order  lag 
plus  dead-time.  These  parameters  are  most  commonly 
obtained  from  a process  reaction  curve  of  the  pro- 
cess. Such  a curve  Is  the  open-loop  response  of 
the  controlled  variable  to  a step-change  In  mani- 
pulated variable.  For  the  reactor,  a first-order 
plus  dead-time  model  was  obtained  by  a method  pro- 
posed by  Miller  [6]  and  the  parameters  obtained 
are  as  follows: 


t f + T(t)  - KWc(t-to) 

K - -0.0333  *F/(lbs/mln) 

t ■ 13.95  min 

t ■ 2.5  min 
o 
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longest  time  constsnt  of  the  process  or  dominant 
pole.  The  derlvetlve  time  Is  considered  as  a 
second  zero  that  compenaates  for  the  second  longest 
process  time  constant.  With  the  two  time  parameters 
thus  determined  to  Insure  fast  loop  response,  the 
gain  can  then  be  adjusted  to  meet  any  specified 
response  criteria.  Gain  correlations  were  obtained 
for  5Z  overshoot  on  the  response  to  a step  change 
In  set-point. 

Martin's  approach  requires  the  first-order  plus 
dead-time  model  parameters  to  tune  a proportional- 
integral  (PI)  controller,  and  the  second-order  plus 
dead  time  parameters  to  tune  a proportional- 
integral-derivative  controller  (P1D) . 


The  process  reaction  curve  was  obtained  by  a step 
reduction  In  cooling  water  rate  of  120  lbs/mln. 
Although  the  parameters  vary  with  the  magnitude  and 
direction  of  the  step,  we  felt  that  it  would  be 
mreallstic  In  an  Industrial  situation  to  perform 
multiple  tests  to  obtain  an  average.  A single  pro- 
cess reaction  curve  Is  difficult  enough  to  get. 

TUNING  METHODS 

Rovira  [1]  and  Lopez  [2]  based  their  tuning  corre- 
lations on  the  first-order  plus  dead-time  model  of 
equation  1.  They  used  parameter  search  techniques 
to  determine  the  controller  parameters  that  result 
in  a minimum  of  the  following  error  functions: 


perature  Controller 


Using  the  parameters  obtained  from  the  process 
reaction  curve  for  the  reactor,  the  PI  controller 
parameters  computed  from  the  various  correlations 
are  shown  in  Table  I.  Note  that  the  parameters 
from  the  Rovira  correlation  are  essentially  the 
same  - within  control  parameter  accuracy  - as  those 
from  the  loop-compensation  approach  for  5Z  over- 
shoot, while  the  Lopez  and  Zlegler-Nlchols  correla- 
tions result  In  tighter  control  parameters.  The 
equation  for  the  PI  controller  Is  given  by: 


m - mQ  + K,  (e  + /edt) 


(1) 


IAE  - rQ  | e | dt 


where  e ■ (set-point)  - T 


ITAE  - /“  |e| tdt 

where  |e|  Is  the  absolute  value  of  the  error  or 
instantaneous  difference  between  the  controlled 
variable  and  the  set-point.  Lopez  also  considered 
the  integral  of  the  squared  error  (ISE)  but  the 
resulting  tuning  parsmeters  gave  highly  oscillatory 
responses. 

While  Lopez  minimized  the  error  Integrals  for  step 
changes  In  disturbance  - assuming  the  process 
dynamics  to  disturbance  were  Identical  to  the  dynam- 
ics to  the  controller  output  signal  -,  Rovira  con- 
sidered step  changes  In  set-point,  a more  demanding 
case.  This  resulted  in  more  conservative  tuning 
parameters. 

Although  Ziegler  and  Nichols  [3]  did  not  base  their 
turning  correlations  on  the  parameters  of  equation  1, 
the  graphically  defined  parameters  from  the  process 
reaction  curve  can  be  converted  to  those  of  the 
first-order  plus  dead-time  model.  They  obtained 
their  correlations  by  empirical  methods  and  based 
them  on  a quarter-decay  ratio.  Unfortunately,  the 
quarter-decay  ratio  Is  a more  oscillatory  response 
than  la  usually  acceptable  In  an  Industrial 
environment. 

These  authors'  loop- covenant  Ion  approach  [4]  con- 
sists of  looking  at  the  controller  as  a dynamic 
compensator  for  the  other  major  covenant  of  the 
loop:  the  process.  In  his  approach  the  Integral 

time  Is  considered  a zero  of  the  controller  trans- 
fer function  that  Is  used  to  compensate  for  the 


The  parameters  for  the  PID  controller  are  given  In 
Table  II.  Note  that  the  Integral  and  derivative 
Limes  for  the  loop-compensation  method  are  of  the 
order  of  magnitude  of  the  process  time  constants. 
Again  the  Lopez  and  Zlegler-Nlchols  correlations 
result  In  higher  gains  and  shorter  Integral  time. 

The  parameters  of  Table  II  are  for  the  following 
PID  controller  equation: 

m - mo  + Kc  (e  + /edt  + T<1  £)  (4) 

COMPARATIVE  SET-POINT  RESPONSES 

The  responses  of  the  PI  controller  to  a 4 *F  rise 
In  set-point  are  given  In  Figure  3.  The  figure 
illustrates  the  closeness  of  the  response  between 
the  Rovira  and  loop-compensation  techniques.  Lopez' 
tuning  la  not  shown  because  It  Is  not  Intended  for 
set-point  changes.  The  quarter-decay  Zlegler- 
Nlchols  response  Is  obviously  too  oscillatory. 

Note  that  the  loop-compensation  response  exceeds 
the  5Z  overshoot  for  which  It  was  designed.  This 
Is  because  of  the  nonlinear  nature  of  the  reactor 
which,  at  this  higher  temperature,  exhibits  a 
higher  gain  than  measured  by  the  process  reaction 
curve.  The  nonlinear  effect  la  evident  In  Figure  4, 
which  shows  the  response  of  the  PI  controller  to  a 
4 *F  drop  In  set-point.  At  the  lower  temperature 
the  reactor  gain  la  lower. 

The  corresponding  responses  of  the  PID  controller 
tuned  by  the  various  correlations  are  shown  In 
Figures  S and  6 for  a 4 *F  rise  and  drop  In  aet- 


point,  respectively.  The  Rovlra  IAE  tuning  h«s  a 
faster  rise-time  than  loop-compensation  tuning.  It 
must  be  kept  in  mind,  however,  that  the  gain  of  the 
controller  can  be  field  adjusted  under  the  loop- 
compensation  concept,  since  it  is  determined  inde- 
pendently of  the  Integral  and  derivative  times. 

RESPONSES  TO  LOAD  INPUT 

The  PI  controller  responses  to  a 20Z  decrease  in 
reactant  feed  rate  are  shown  in  Figure  7.  In  this 
case  the  Lopez  tuning  is  Included  since  it  is 
Intended  for  load  changea.  Rovlra' s tuning  is 
essentially  the  same  as  loop-compensation  tuning 
for  the  PI  controller  as  mentioned  earlier.  Note 
the  oscillatory  behavior  of  the  Zlegler-Nlchola  and 
Lopez  responses.  Again  it  appears  that  the  loop- 
compensation  gain  could  be  adjusted  upwards  to 
speed  up  its  response.  The  corresponding  PID  con- 
troller responses  are  shown  in  Figure  8.  The  Lopez 
and  Zlegler-Nlchols  responses  are  less  oscillatory 
in  this  case.  Again,  loop-compensation  is  rather 
conservative,  and  seeme  that  it  could  stand  a 
higher  gain. 

SUMMARY  AND  CONCLUSIONS 

Tuning  techniques  based  on  simple  models  of  the 
process  have  been  shown  to  result  in  very  satis- 
factory responses  when  applied  to  the  temperature 
control  of  a nonlinear  chemical  reactor.  The 
tuning  methods  of  Rovlra  and  loop-compensation  pro- 
duced the  best  responses  to  set-point  changes  in 
terms  of  stability  and  speed  of  response.  In  gen- 
eral, set-point  changes  are  more  demanding  in  con- 
troller performance  than  load  changes,  and  are 
equivalent  to  load  changes  when  the  process  dynam- 
ics to  disturbance  inputs  are  much  faster  than  its 
dynamics  to  the  manipulated  variable. 

Although  conservative  in  terms  of  responses  to  load 
changes,  the  fact  that  the  loop-compensation  con- 
cept decouples  the  adjustment  of  the  gain  from  that 
of  the  time  parameters  allows  the  operator  to 
obtain  any  desired  response  criteria  by  adjustment 
of  only  one  knob.  This  is  not  the  case  for  the 
other  methods  for  which  all  of  the  controller  param- 
eters are  Interrelated. 
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NOTATION 

b Second-order  model  damping  parameter 

c Second-order  model  frequency  parameter 

e Controller  arror 

IAE  Integral  of  the  absolute  value  of  the  error 
1TAE  Integral  of  time  averaged  abaolute  value  of 
Che  error 
K Model  gain 

K,.  Controller  gain 

m Controller  output  elgnal 

■o  Initial  controller  output  elgnal 

t Tima 

t„  Modal  dead-time 


T Temperature  of  the  reacting  fluid 

Tg  Controller  derivative  time 

Tjl  Controller  Integral  time 

W Mass  rate  of  feed  (load  input) 

Wc  Maas  rate  of  cooling  water  (manipulated 
variable) 

t Model  time  constant 
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TABLE  I 

PI  CONTROLLER  TUNING  PARAMETERS  - BACKMIX  REACTOR 


Techniques 

*c 

(lbs/mln/'F) 

tun 

Loop-Compensstlon 
(5Z  overshoot) 

-0.584 

13.95 

Rovlra  (Set-Point) 
(IAE) 

-0.666 

14.50 

Lopez  (Load) 
(IAE) 

-1.072 

6.805 

Zlegler-Nlchols 

-1.003 

8.32 

(quarter  decay) 


Zlegler-Nlchols 


Zlegler-Nlchols 


Figure  7.  Closed- loop  response  of  FI  controller  for  s 
20Z  decrease  In  reactant  feed  rate 


Closed-loop  response  of  FID  controller 
for  a 4°F  drop  In  set-point 


Figure  6 


Nichols 


Loop-Compensation 


20  30  40 
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Figure  8.  Closed-loop  response  of  PID  controller 

for  a 202  decrease  In  reactant  feed  rate 
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